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Functions which are covariant or invariant under the transformations of a compact 
linear group G acting in a euclidean space IR", can be profitably studied as func- 
tions defined in the orbit space of the group. The orbit space is the union of a 
finite set of strata, which are semialgebraic manifolds formed by the G-orbits with 
the same symmetry. In this paper we provide a simple recipe to obtain rational 
parametrizations of the strata. Our results can be easily exploited, in many phys- 
ical contexts where the study of covariant or invariant functions is important, for 
instance in the determination of patterns of spontaneous symmetry breaking, in 
the analysis of phase spaces and structural phase transitions (Landau's theory), 
in covariant bifurcation theory, in crystal field theory and in most areas of solid 
state theory where use is made of symmetry adapted functions. An example of 
utilization of the recipe is also discussed at the end of the paper. 



1. Introduction 

The determination of properties of functions which are covariant or invariant 
under the transformations of a compact linear group (hereafter abbreviated 
in CLG) G is often a basic problem to solve in many physical contexts. G- 
invariant functions play an important role, for instance, in the determination 
of patterns of spontaneous symmetry breaking and structural phase transitions 
(Landau's theory), in covariant bifurcation theory, in crystal field theory and 
in most areas of solid state theory where use is made of symmetry adapted 
functions. 

It will not be essentially restrictive, in the following, to assume that G is a 
matrix subgroup of 0„ (M) ; it will be thought to act in M" , thus defining the 
G-space (G,IR"). Complex linear groups can always be transformed into real 
linear groups through a process of realification and compact real linear groups 
are equivalent to orthogonal linear groups. 

A possible approach to the study of invariant functions, that fully exploits 
invariance and regularity properties ^, is based on the fact that a G- invariant 
function f{x), x G IR", is a constant along each G-orbit and can, therefore, be 
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considered as a function / in the orbit space of (IR"/G) ^. 

Geometric invariant theory ^ suggests how to get, in principle, / and 1R"/G. 
In fact, theorems by Hilbert, Schwarz and Nagata ^ assure that a polynomial or 
C°° — , G-invariant function /(.x), x E IR", can be expressed as a polynomial or, 
respectively, C°° function /(p) of a finite minimal set p{x) = {pi{x), . . . ,Pn{x)) 
of polynomial G-invariant functions {minimal integrity basis of the ring of G- 
invariant polynomial functions (hereafter abbreviated in MIB)). The elements 
of a MIB separate the orbits of G. 

Depending on G, there may be relations among the elements of a MIB 
{non-coregular G): 

= 0, Va; e ]R" and A= 1,..., . (1) 

The associated equations 

Fa{p)=0, A=1,...,, (2) 

define an algebraic variety in IR"^, that is called the variety 2 of the relations 
(among the elements of the MIB). 

The range p{Si^) of the orbit map p provides a diffeomorphic image of K"/ G 
in Z C IR'?. The set p{iSC) is a connected semi- algebraic subset (i.e., a subset 
defined by algebraic equalities and inequalities) of IR'', formed by the union 
of semialgebraic (connected) manifolds of different dimensions, called primary 
strata. The G-orbits with the same isotropy type form semialgebraic sets 
(isotropy type strata), whose connected components arc primary strata. Almost 
all the points of IR"/G belong to a unique stratum, the principal stratum, which 
is a connected open dense subset of 'iBT' /G. The boundary of the principal 
stratum is the union of disjoint singular isotropy type strata. All the strata 
lying on the boimdary of an isotropy type stratum are open in its topological 
closure. The following partial ordering can be introduced in the set of all the 
orbit types: {H\ < [K], if H is conjugated with a proper subgroup of K. The 
orbit type of an isotropy type stratum is contained in the orbit types of all 
the strata lying in its boundary; therefore, more peripheral strata are formed 
by orbits with higher symmetry under G transformations. The number of 
distinct orbit types of G is finite and there is a unique minimum orbit type, 
the principal orbit type, pertaining to the principal stratum; there is also a 
unique maximum orbit type [G] , corresponding to the image through tt of all 
the points of H", which are invariant under G. 



^The orbit space R"/G is defined as the quotient space obtained through the equivalence 
relation between points belonging to the same orbit), endowed with the quotient topology 
and diflferentiable structure. 
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The principal stratum is connected, open and dense in ]R"/G. The orbit 
type [H] of a stratum^ is contained in the orbit types of all the strata lying in 
its boundary; therefore, more peripheral strata of IR"/G are formed by orbits 
with higher symmetry under G, the principal orbit type is minimum and the 
maximum orbit type is [G] , corresponding to the image through the canonical 
projection, tt : H" ]R"/G, of the set of points of K", which are invariant 
under G. This set contains at least the origin of M". 

The following theorem offers a simple recipe to derive the algebraic equa- 
tions and inequalities defining the range of an orbit map and its strata ^ (see 
also 

Let P{p), p e IR' denote the q x q real symmetric matrix defined through 
the relations 

n 

Pab{p{x)) = J2 ^jP»i^) djPbix), \/x e R". (3) 

Then, p(IR") — {p ^ Z \ P{p) > 0} and the k-dimensional primary strata 
are the connected components of the semialgebraic set {p ^ Z \ P{p) > 
0, rank(P(p)) = k}. 

If G is coregular, the Pi's, with range in the semialgebraic set A defined by 
the inequalities P{p) > (and no equation!) provide a one-to-one parametriza- 
tion of the principal stratum of ]R"/G. If the stratum is singular, the pi's are 
not independent parameters, being bounded also by the equations defining the 
stratum. In applications or, simply, to get a better understanding of the ge- 
ometry of the stratum (its connectivity properties, its boundary, etc.), one 
often needs to solve explicitly these equations, which are obtained in implicit 
form from the proposition recalled above. Often, this cannot be done using 
standard algorithms. 

In this paper we propose a general method to derive the equations defining 
the singular strata of the orbit space of any CLG in the form of explicit rational 
parametric relations. In the next Section we shall present and prove our results 
and in Section 3, we shall illustrate them, by working out an example. 

2. Parametrization of strata in orbit spaces of compact linear 
groups 

In this section, we shall prove that the possibility of parametrizing the prin- 
cipal stratum of the orbit space of a coregular group in terms of parameters 



^[H] is the class of subgroups of G conjugate with H and, by definition, is the class of 
isotropy subgroups of G at points of the stratum. 
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related to the elements of a MIB can be extended to singular strata. The re- 
sult stems from the proof of the following statement: The set of internal points 
of the topological closure of an isotropy type stratum, with orbit type [H], is 
diffeomorphic to the principal stratum of the orbit space of a group space of the 
stabilizer of H (in G). If the stabilizer turns out to be coregular, the construc- 
tive proof of this proposition will provide a one-to-one rational parametrization 
of the stratum. 

Since, in view of possible applications, we are mainly interested in the possi- 
bility of getting a one-to-one parametrization of a stratum, in the presentation 
of our results, stress will be laid on this aspect of the matter. 

Let G and {p} be defined as in the previous section (no coregularity as- 
sumption will be done on G) and if be a proper isotropy subgroup of the 
G-space (G,]R"). We shall denote by S the stratum of (G,]R") with orbit 
type [H] and by S its image in the orbit space IR"/G. In the following, we 
shall identify IR"/G and its strata S with their images p(IR") and p{S). Let 
us, also, define 

y = {a; e ]R" I G^ D H}, 

(4) 

v = {xemj'\G^ = H}. 

The non trivial sets V and V have the following properties, which are more 
or less immediate consequences of their definitions: 

(1) V = {x G H" \ hx = x,h G H} is a linear subspace of IR", let us call u 
its dimensions; 

(2) V is the topological closure of V: V = V; 

(3) V = S (IV, so that every G-orbit lying in S has at least a point in V 
and S = G - V, S = G -V, where the bar denotes topological closure, 
so that 

p{V) = S, p{V) = S (5) 

and 

rank(P(ar)|^^^) = dim(5), rank(P(a;)|^^^^^) < dim(5). (6) 

A G-orbit of S may intersect V in one or more points. In the first case, every 
G-orbit of S intersects V in a point and. owing to item 3, the coordinates of the 
points of V provide a one-to-one rational (in effect, linear) parametrization of 
the orbits of G lying in S, obtained by solving the system of linear equations 
hx = X yh G H. The allowed range of x, x G V, is determined by the 
inequalities assuring that rank(P(a;)) =dim(S'). The parameters x, with range 
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V, would provide, in this case, the one-to-one parametrization of the stratum, 
thus solving our problem. 

In general, however, the intersection of a G-orbit, of orbit type [H], with 
the stratum 5* docs not reduce to a single point. So, a sounder analisys is 
required, that will be the object of the rest of this section. 

Two distinct points, x and gx, g G G, of the same G-orbit, lay in V iff 
hgx = gxVh E H. This is equivalent to the condition that g is an element of 
the stabilizer Stab(ir, G) of H in G: 

Stah{H, G) = {gGG\ gHg-^ = H}. (7) 

The intersection of a G-orbit of S with V is, therefore, the Stab(iJ, G)-orbit 
through X. 

This result can be refined. In fact, it is easy to realize that the relation 

hgx = gx is satisfied for all x E V and all h E H, iS g E Stab(iJ, G). 
Therefore, Stab(ff, G) is the largest subgroup of G under which V is stable; a 
known result. 

In the group space (Stab(fl', G),V), the isotropy subgroup at a point of gen- 
eral position is H (or, more precisely, the linear group obtained by restriction 
to V of the action of H). Therefore, the principal isotropy type is [H] = {H} 
and the principal stratum S satisfies the following relation: 

1^ D S D V. (8) 

It has to be noted that, the principal stratum S could contain V in a strict 
sense, since, at points of V\V, the conjugacy class of the isotropy subgroup of 
Stab(i?, G) could be smaller than the conjugacy class of the isotropy subgroup 
of G. 

H is an invariant subgroup of Stab(i?, G) and it is a subgroup of all the 

isotropy subgroups of (Stab(iJ, G), y). Thus, the restriction of the action 
of Stab(if, G) to V defines a linear group K (and the group space {K,V)), 
isomorphic to the quotient group Stab(i/, G)/H, through the relation 

{sH)v^sv, s EStah{H,G), sH EStah{H,G)/H. (9) 

So, we can conclude that the intersection of a G-orbit of S with V is an orbit 
of {K, V). Obviously, the orbit spaces V/K and V/Stah{H, G) are isomorphic 
and can be identified. 

We can rephrase the result just obtained, by claiming that the points of 
S are in a one-to-one correspondence with the points of a, possibly proper, 
subset V of the principal stratum S of the orbit space V/ K and the following 
relations hold true: 

^ = P{V) D p(S) D p(V) = S. (10) 
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Since, as stressed, we know how to parametrize a principal stratum, our prob- 
lem is reduced to the determination of V. We shall show that this can be easily 
done, making use of (6). 

To make simpler the solution of the problem, let us assume that an or- 
thonormal basis has been introduced in H" such that the first v elements of 
the basis yield a basis for the vector space V . Then, if we denote by Vj_ the 
orthogonal complement of V in IR", the subspacc Vj_ is stable under [K and) 
iJ, owing to the orthogonality of the transformations of G. Since V contains 
all the if-invariant vectors of M" , there is no non trivial if-invariant vector in 

To attain our goal, we shall start from a convenient parametrization of 
the principal statum S oiVjK in terms of I real parameters A, related to a 
MIB (Ai(u), . . . , A;(t;)) of the ring of polynomial /T-invariant functions. This 
parametrization will be global and one-to-one, if K turns out to be corcgular. 
In this case, the range of A has to be restricted to the positivity region of the 
P-matrix A(A), associated to the MIB {A}: 

V 

Aa/3(A(t;)) = Aa^(w) = ^ di\o,{xx, • • • , a;^) diX/3{xi, x^). (11) 

i=l 

As recalled in the previous Section, in fact, the orbit space V/K and its prin- 
cipal stratum S can be identified, respectively, by the semialgebraic sets X{V) 
and A(S): 

g = A(S) = {A e I A(A) > 0}, S = A(y) = {A e ]R' I A(A) > 0},(12) 

the second set being the closure of the first one. Moreover, the definition of V 
and equation (8) imply 

V = A(V) C A(S) = S. (13) 

If K is not coregular, only a local one-to-one parametrizations can be 
obtained for S, by eliminating redundant elements in the set of parameters 
(Ai,...,A(), through the solution of the algebraic relation(s) among the el- 
ements of the MIB {A}, and imposing convenient semi-positivity and rank 
conditions on the matrix A(A). 

A one-to-one local or global parametrization of S yields, obviously, also 
a local or global one-to-one parametrization of V, provided that additional 
restrictions are imposed on the range of A, whenever V is a proper subset of 
E. The correct bounds can be obtainc^d in the following way. 

We shall only consider the case of a coregular K, the extension of the results 
to non coregular K's will be straightforward, but, as just stressed, may lead 
to a loss of globality. 
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When X spans V , the elements of the MIB {pi{x), . . . ,pq{x)) of G define 
a set of ii'-invariant polynomial functions of v = {x\, . . . ,Xv) € V. Therefore, 
by the Hilbcrt theorem recalled in the previous Section, they can be expressed 
as polynomial functions of A, that is, 

p|^ = <^oA. (14) 

Let us remark that, in our assumptions, there are no relations among the 
elements of the MIB {A} and, consequently, possible relations Fa{p) = among 
the Pi's are identically satisfied for p = 0(A). 

^Prom (14) and (8), one immediately obtains 

^ = p{V) = cj>{V/K) D <^(S) D </.(V) = p{V) = S. (15) 

and, since S is the set of internal points of V/K, ^{T,) will coincide with the 
set of internal points of the closure S of S. This set does not coincide with S 
if S contains, in its interior, points representing bordering strata of S. This 
happens if the set S is not connected or not multiply connected. 

The identification of points A e S, whose image 0(A) ^ S can be obtained 
in the following way, using (6). 

Let X = v(Bv± yield the decomposition of a; G IR" in its vector components 
V G V and v± G V±: v = (xi, . . . , x^,), v± = (x^+i, . . . , a;„). Then, a G- 
invariant polynomial f{x), can be thought of as a polynomial in v and v± and 
it is easy to realize that it cannot contain linear terms in v± , being v invariant 
under H. As a consequence, 

dif{x) = 0, for X E V and i = i/ + 1, . . . ,n. (16) 

Now, starting from the very definition of P{x) (see (3)), for every x = v € V 
we obtain, using (16) and the identity p{v) = 0(A(w)), 

Pab{p{v)) = ELi diPaiv) diPb{v) 

= (j{XfA{X)J{X) ) , 

where, the apex T denotes transposition and J(A) is the Jacobian matrix of 
the transformation p = 4>{X): 

Jaa{X) = da(f)a{X), a = l,...,q, a=l,...,l. (18) 

So, we can conclude that, for all A e V/K, 

P(<A(A)) = J(A)TA(A)J(A). (19) 



EiLi diXa{v) diX0{v) 



(17) 
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Equation (19) leads to an easy identification of the points A e V, that is, 
of the points A G S whose image p = 4>{X) is in S. In fact, from (6), these 
points are characterized by the conditions A e S and Ta.nk{P{^{X))) = I. Since, 
for A G S, the matrix A(A) is positive definite with rank I, we can conclude 
that the range of A assuring the location of the point p = (f){X) in S is the 
semialgebraic subset A of 1R', determined by the following inequalities: 

A = {A e I A(A) > and rank(J(A)) = I}. (20) 

For A G A, the relation p = 0(A) yields a global rational parametrization for 
S. 

It will be worthwhile to stress that the parametrization we have siiggcsted 
is, in some way, canonical: The unique arbitrariness in the choice of the pa- 
rameters is related to the choice of the MIB's. 

An important byproduct of the result just proved is a simple test of the 
connection of 5, that could be difficult or impossible to read directly from the 
equations of the stratum in implicit form. In fact, the condition 

boundary of 5 = boundary of S 

is equivalent to rank(J(A)) = I for all A G E. 

If K is not coregular, there are polynomial relations -Fa (A) = among the 
elements of the MIB (Ai, . . . , A;) and / >dim(S'). As already stressed, to obtain 
a one-to-one parametrization of the points of S by means of the A^'s, one has 
to eliminate the redundant parameters by solving the equations -Fa (A) = 0. 
This may be feasible, but only locally and the resulting parametrization will 
not be global and possibly not rational. 

3. An Example 

In this Section we shall show how the parametrization technique works in 
a simple example. The notations will be the same defined in the previous 
sections. 

We shall consider the orthogonal linear group G defined by the action 
of the real group 0(3) in the real eighth-dimensional space spanned by the 
independent components xi, . . . ,X5 of a symmetric and traceless tensor Q and 
the three components X6,xt,xs of a polar vector P. To be specific, if O is a 
generic 3x3 real orthogonal matrix, the transformation rules of the 
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obtained from the following relations; 




/ 



(21) 




a,/3= 1,2,3. 



(22) 



It will be worthwhile to remark that the reflection I = diag(— 1, —1,-1) G 0(3) 
reverses only the signs of the last three coordinates (xg, .tj, xs) and, due to the 
symmetry of the tensor Q, each G-orbit contains points whc^re the tensor Q 
takes on a diagonal form (x^ = X4 = x^ = 0). This remark makes much easier 
the determination of "typical points" in strata, from which the orbit type of 
the stratum can be identified. 

The real orthogonal linear group G, just defined, is coregular (see ^) and 
its isotropy subgroup at a generic point of is trivial. As a consequence, 
the principal orbits are 3-dimensional manifolds, the principal stratum has 
dimensions 5 and there will be five independent elements in a MIB (see, for 
instance ^), which can be chosen in the following way: 
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i=6 



6V2 TtQ^ 

-2V3xi^ + 6 VSxi - 3 V3 .Ti xs'^ -9x2 x^'^ 
—3 ^f^xi Xi' +9x2Xi^ + 18 xa a;4 + 6 V^xi x^^ 



(23) 



3\/2 ^P„Qo;3P/3 



— 2 VSxi X6^ + 6X3X6X7 + VSxi X-j^ — 3X2 X'j^ 
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+6X4,X6 Xa +6X5 X7 Xg + V^Xi Xs"^ + 3X2 Xs"^ , 

= Axi'^xe'^ + 3x3'^ xq'^ +3a;4^a;6^ - 2^3x1 X3 xq xr 

—6 X2 'C;! -'I'd xj + 6 X4 X5 xg X7 + xi^ X'j^ — 2 \f2>x\ Xi Xt^ 
+3 X2^ xt^ + 3 X'i' X'j^ + 3 x^ xi^ — 2 \f2)X\ X4 xe xg 
+6 X2 X4 xg xs + 6 X3 X5 xg xs + 6 X3 X4 X7 xg + 4 -s/S xi X5 X7 xg 
+xi^ xg^ + 2 \/3 xi X2 xg^ + 3 X2^ xg^ + 3 X4^ xg^ + 3 x^ x^ . 

The corresponding P-matrix elements can be easily calculated from their 
very definition (3) 



Pi a 


= 2daPa , 1 < a < 5 


P22 


= 4p2 


P23 


= 


P24 


= 4j34 


P25 


= 4p5 


P33 


= 108 (pi - P2f 


P34 


= 18 {-2pip2 + 2p2^ +P5) 


P35 


= 12 (p2 P3 + Pi Pi - P2 P4) 


P44 


= 12 (pa'+Ps) 


P45 


= 4 {P2P3 + 3piP4- P2P4) 


P55 


= 1 {P3P4 +P4' + QpiPs) , 



where the set {^1,^2,^3,^4,^5} = {2,2,3,3,4} specifies the degrees of the 
polynomials of the MIB. 

The determinant of the matrix P{p) factorizes and only one of the two real 
irreducible factors, that we shall call A{p), turns out to bo active ®, that is, to 
be relevant in the determination of the boundary of the orbit space IR^/G: 

A{p) = 3p2^P3^ + 18piP2^P3P4 - 18p2'^P3P4 + 27pi'^p2Pi^ 

-MpiP2^P4^ + 27p2^P4^ +P3PA^ - 9P2P3P4P5 (25) 
-9piP4^P5 + 9p2P4^P5 - 27piP2P5^ + 27^2^^52 + 9^5^ • 

The relations assuring that P{p) > and has rank 4, define a unique 4- 
dimensional stratum S^'^^ in the orbit space. Using well known matrix theory 
results, these conditions could be explicitly written, for instance, in the form 
A{p) = and Mi{p) > 0, i = 1, ... ,4, where Mj is the sum of the order i 
principal minors of the matrix P{p): A frightening set of conditions that it is 
not worthwhile to write down explicitly. 
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Even in this simple example one immediately realizes that the choice of a 
minimal set of explicit algebraic relations providing a cylindrical decomposition 
for the semi-algebraic subset S*^^^ of IR^ would be a really hard task (for the 
more peripheral strata, instead, the problem is much easier to solve). An 
immediate application of the results proved in the previous Section, on the 
contrary, leads to a simple rational global parametrization of the sub-principal 
stratum 

A "typical point" of the stratum is Xt = (1,1,0,0,0,0,1,1). The corre- 
sponding isotropy subgroup of G is the TL^ group generated by the reflection 

diag(l,l,-l,-l,l,-l,l,l)eG, 

representing the element diag(— 1, 1, 1) e 0(3). 

The vector space V formed by the if-invariant vectors of H" turns out to 
be be the following: 

The elements of 0(3), corresponding to elements of Stab(iy, G), are 
block-diagonal matrices of the form diag(±l,0), with O € 0(2). There- 
fore, Stab (7J, G) is isomorphic to a group TLi®0(^ and the quotient group 
K = Stab {H, G) /H is the representation in (induced by the representa- 
tion G of 0(3)) of the Oi(2) subgroup of 0(3), formed by the rotations around 
the first axis. 

Using coordinate v = {xi,X2, X5,X7, xs) for a vector v gV, the elements of 
K turn out to be block diagonal matrices leaving stable the following subspaces 
V^,1<A<3: 

= {(xi, 0,0, 0,0) e V\xi e IR}, 

1/2 = {(0, X2, X5, 0, 0) e V\X2 , Xg e ]R} , 

= {(0, 0, 0, X7, Xs) G V\X7, , xg € ]R} . (26) 

A proper rotation r{(j)) G Oi(2) of an angle (j) and the reflection diag(- 
l,l)e Oi(2) are represented in the primed basis of V by the 5x5 matrices 
h{(j)) = 1 ® r(-20) © r(0) and j = diag(l, -1, 1, -1, 1). 

Noting that the complex variables Zi = X2+ix^ and Z2 = Xy-t-ixg transform 
into e~2*'^zi and, respectively, e'*'z2 under a rotation h{<f)) and into zl and, 
respectively, —Z2 under a reflection j, it is easy to realize that a possible MIB. 
for {H, V) is the following: 



12 



Ai = xi , 

A2 = = xl+xl, 

A3 = |22p = Xt+xI, 

A4 = 2Re(2;i^;2^) = 2 (a;2a;7^ - 2a;5a;7a;8 - X2Xs^) ■ (27) 

It is, now, easy to express the p's in term of the A's, p = ^(A), and to 
check that the following expressions, obtained in this way, identically satisfy 
the equation A{p{X)) = (see (25)) entering the definition of the stratum 5^^^: 



0i(A) = A12 + A2 + A3, 

02(A) = A3 , 

03(A) =2 V3 {Xi^- 3X1X2) , 

04(A) = (2AiA3 + y3A4) , 

05(A) = Ai2A3 + 3A2A3 + VsXiXi . 



(28) 



As explained in the previous section, since the group K is coregular, the 
range A for A, is the region where the P-matrix, A(A) associated to the MIB 
{A} is positive definite and the rank of the jacobian matrix J(A) of the trans- 
formation 0(A) is maximum (=4). The matrices A(A) and J(A) are easily 
calculated to be 



A(A) = 



/I 

04A2 2A4 

4A3 4A4 
V02A4 4A44 (A3^ + 4A2A3)/ 



(29) 



and 



/ 2Ai 1 1 \ 

10 

6 a/3 (Ai^ - A2) -6\/3Ai 
-VSXs -VSXi -f 

V 2AiA3 + a/3A4 3A3 Ai2 + 3A2V3Aiy 



(30) 



The conditions assuring the positivity of A(A) > can be written in the form 

A = {(Ai, A2, A3, A4) e I A2 > and A3 > and A^ < 4A2A3^} (31) 

and the rank of J(A) turns out to equal four for all A e A. So we can conclude 
that the stratum is connected, a piece of information that would be difficult to 
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derive directly from the relations defining the stratum in implicit form {{A(j}) = 
and Mi{p) > 0). 

The paramctrization obtained for the sub-peripheral stratum S'^^-') turns 
out to be useful, also, because bordering values of A immediately determine 
bordering values of p, corresponding to more peripheral strata of IR^/G. 
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